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Neutron powder diffraction has been used to investigate the structural and magnetic behavior of 
the isoelectronically doped Fe pnictide material PrFei-xRu^AsO. Substitution of Ru for Fe sup- 
presses the structural and magnetic phase transitions that occur in the undoped compound PrFeAsO. 
Contrary to the behavior usually observed in 1111 pnictide materials, the suppression of both the 
structural and magnetic transitions does not result in the emergence of superconductivity or any 
other new ground state. Interestingly, PrFeAsO itself shows an unusual negative thermal expansion 
(NTE) along the c-axis, from 60K down to at least 4K; this does not occur in superconducting 
samples such as those formed by doping with fluorine on the oxygen site. We find that NTE is 
present for all concentrations of PrFei_ x Ru x AsO with x ranging from 0.05 to 0.75. These results 
suggest that the absence of superconductivity in these materials could be related to the presence of 
NTE. 



I. INTRODUCTION 

Since the discovery of Fe-based superconductors in 
2008 |TJ [2] there has been enormous interest in obtain- 
ing a detailed understanding of the properties of closely 
related materials. Studies of substitutionally altered or 
doped materials have proved particularly valuable and 
have led to the discovery of several superconducting fam- 
ilies. Here we report new results for the isoelectronically 
doped 1111 pnictide material PrFei-^Ru^AsO. 

The 1111 FeAs based materials (typified by LaFeAsO) 
share several common characteristics with most of the 
other Fe-based superconductors. As the temperature is 
lowered the undoped parent material does not reach a su- 
perconducting ground state, but undergoes a structural 
transition, followed or accompanied by an antiferromag- 
netic spin density wave (SDW) transition of the moments 
associated with the Fe atoms and sometimes antiferro- 
magnetic (AFM) ordering of the rare earth moments [2]. 
For superconductivity to emerge in the material, all of 
these transitions must be eliminated or suppressed down 
to a sufficiently low temperature [3[ [4]. This can be done 
by doping with the appropriate element, e.g. F on the O- 
site in CeFeAsO[5], Co on the Fe-site in LaFeAsO [6j 0, 
etc. In certain members of the Fe oxypnictides, supercon- 
ductivity can also be induced by simply applying pres- 
sure externally, e.g. LaFeAsO under 120kbar for a Tc of 
21K[8] and BaFe 2 As 2 under 40kbar for a T c of 29K0. 
The application of pressure does not always produce this 
result, for example CeFeAsO achieved no superconduc- 
tivity under pressures up to 500kbar[T0]. 

There are situations where AFM ordering and super- 
conductivity can coexist, for example in 122 materials 
such as in BaFe2- ;E Co :E As2[II]. Although this occurs 
in several 122 materials [12] the general presumption 
that superconductivity and AFM ordering are competing 



ground states remains valid [T3 ] H4 ] . Notably coexistence 
is much less common in 1111 materials such as LaFeAsO, 
where long range AFM ordering is usually completely de- 
stroyed before superconductivity can emerge, although 
there are some exceptions such as SmFeAsO^Fi - x [T5] . 
The suppression of the structural and magnetic transi- 
tions can also lead to new non-superconducting ground 
states, for example ferromagnetism[T6] . It has been 
concluded that suppressing the structural and magnetic 
transitions is usually necessary, but not sufficient, for in- 
ducing superconductivity [3]. 

The unit cell of the 1111 parent compound PrFeAsO 
is illustrated in Figure 1. Most of the materials with 
the highest reported superconducting T^'s are doped 
on either site of the rare-earth-oxygen layer [2] [14]. Su- 
perconducting Tc's are generally lower for materials 
doped on the Fe-As layer, possibly due to undesir- 
able lattice distortions and disorder in the conducting 
Fe-As layer [T7] [18]. To better understand this, iso- 
electronic substitution on the Fe site has proved very 
instructive since unlike carrier doping it enables the 
study of materials without directly altering the Fermi 
level [T9]. Moreover, isoelectronic doping with Ru can 
lead to superconductivity in certain Fe-oxypnictides, 
such as the 122 materials, BaFeo-^Ru^ As9 [T7] and 
SrFe9_^Ru^ Asq [20] . In contrast, however, Ru/Fe sub- 
stitution in 1111 materials has not been observed to lead 
to superconductivity even after the complete suppression 
of the structural and magnetic transitions: the studied 
systems include La(Fe,Ru)AsO[4], Pr(Fe,Ru)AsO[l9], 
Sm(Fe,Ru)AsO and Gd(Fe,Ru)AsO|2T]. 

McGuire et al. showed that Ru doping in PrFeAsO 
does not induce superconductivity down to at least T = 
2K[19]. Even when the structural and magnetic transi- 
tions were suppressed completely there was no indication 
of any new ground state [19]. Bulk and transport mea- 
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FIG. 1. The PrFei-xRu^AsO room temperature tetragonal 
unit cell showing the rare-earth oxygen (top, bottom) and 
Fe-As (middle) layers, the Fe-As pyramid and the Fe-As-Fe 
angle <j>. Note that in the orthorhombic structure there are 
two inequivalent values for <fi. 
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surements on PrFe i-^Ru^ AsO showed hints of anomalies 
associated with a possible phase transition for concentra- 
tions up to x = 0.5 or 0.6[19]. This is a rather high level 
when compared with carrier doped Fe-pnictides, where 
the transitions are suppressed with less than 10% doping 
in most cases. Thus, Ru-doping provides a wide stoichio- 
metric window for studying the progressive suppression 
of the structural and magnetic transitions. Ru-doping 
studies have been also carried out on F-doped supercon- 
ducting samples, such as LaFei_^Ru^AsOn aoFn n [22] 
and SmFe1_aRuaAsO0.85F0.15 [23]. The general results 
show that the superconductivity in these samples is not 
sensitive to isoelectronic doping except at very high levels 
of the order of x = 0.3 and above. 

A phenomenon of particular interest in the parent com- 
pound PrFeAsO is the observation of negative thermal 
expansion (NTE) along the c-axis, where the sample 
gradually expands along the c direction as it is cooled 
below 60K[24 . Notably such NTE is absent in super- 
conducting F-doped PrFeAsO [24 . Recent results have 
shown that NTE also occurs in NpFeAsO were it is as- 
sociated with magnetic ordering of the Np ions [25]. A 
different NTE behavior has been observed in the su- 
perconducting samples of Co doped BaFe2As2 where a 
sudden onset of NTE along the c axis appears at Tc- 
This feature is absent in the non-superconducting parent 
BaFe2As2[26 . It seems worthwhile to explore whether 
or not there is a systematic correlation between the NTE 
and the presence or absence of superconductivity. De- 
spite voluminous reports of the dependence of lattice con- 
stants on stoichiometry in Fe-As based materials there 
has been surprisingly little published data showing the 



FIG. 2. (a): Heat capacity vs T for PrFe i-^Ru^ AsO, the 
peak at 14K indicates the AFM ordering of Pr in PrFeAsO. 
No peak is observed for x > 0.1; (b):x _1 vs T, the inset shows 
a closer look at the lower temperature range for x vs T. 



temperature dependence of the lattice constants, espe- 
cially the c-axis. 

In this paper we report new results on Ru doped 
PrFei-^Ru^AsO, , including detailed studies of the tem- 
perature dependence of lattice parameters. As mentioned 
above undoped PrFeAsO shows NTE. As discussed be- 
low, we observe that Ru-doping in PrFeAsO suppresses 
the structural and magnetic phase transitions without 
leading to superconductivity; however NTE along the c 
axis persists for all doping levels up to at least 75% Ru on 
the Fe site. This paper is organized as follows: Sample 
synthesis and preliminary characterization is described in 
section II, the results of neutron scattering experiments 
are described and discussed in section III, with conclu- 
sions in section IV. 



II. SAMPLE SYNTHESIS AND 
CHARACTERIZATION 

Polycrystalline samples of PrFe i-^Ru^ AsO with x = 
0.05, 0.1, 0.33, 0.4, 0.5, 0.6 and 0.75 were synthesized us- 
ing methods reported previously by McGuire et al [19]. 
In each case, a stoichiometric mixture of powders of 
PrAs, Fe203, Ru02, Fe and Ru was crushed and mixed 
thoroughly in a glove box, pressed into pellets and re- 
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FIG. 3. (Upper): Fully refined neutron powder diffraction 
pattern collected at HB2A of HFIR for PrFei_ x Ru x AsO at 
4K; (Lower): 220jv peak for x = 0.1, 0.33 and 0.75 at T = 
200K (red circles) and 4K (black squares). Peak splitting is 
clearly visible only for x=0.1. 



ceived two firings at 1 200-1 250° C. Transport property 
data and results of electronic density of state calcula- 
tions for the same series of materials have been reported 
previously [19]. The undoped parent compound PrFeAsO 
undergoes a tetragonal-orthorhombic structural transi- 
tion at T ~ 150 K, followed by a SDW transition at 
T ~ 140 K associated with the AFM ordering of Fe, 
and then finally AFM ordering of the Pr moments at 
14 K[27]. Previously published transport measurements 
for PrFei-^Ru^AsO have shown evidence for the sup- 
pression of long range order of the Pr moments for x > 
0.1 and structural/SDW transitions for x > 0.5[19 . 

Prior to the neutron powder diffraction experiments, 
we have carried out heat capacity and magnetization 
measurements. As seen in Fig. 2(a) and (b), the AFM 
ordering of Pr at 14K in the x=0 material is indicated 
by a peak in the heat capacity and a kink in the magne- 
tization. These features disappear for x > 0.1. 



III. EXPERIMENTAL RESULTS AND 
DISCUSSION 

We performed neutron powder diffraction on 
PrFei-^Ru^AsO using the neutron powder diffrac- 
tometer HB-2A (for x = 0.1, 0.33 and 0.75) and the fixed 
incident-energy triple- axis HB-1A (for x = 0.05 and 0.1) 
of HFIR at ORNL to study the structural and magnetic 
transitions respectively. Additional measurements were 
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FIG. 4. The ratio Xtetra given by the Rietveld refinement 
with the P4/nmm (tetragonal, space group 129) structure, 
to xlrtho obtained from refining to the Cmma (orthorhombic, 
space group 67) structure. Any value larger than 1.0 indicates 
the Cmma model provides a better fit than P4/nmm. The 
represents data collected from both HB-2A and POWGEN. 



performed at the high-resolution powder diffractometer 
POWGEN of SNS at ORNL to study the effect of Ru 
doping on the anomalous NTE and also to have a closer 
look at the suppression of the structural transition (for 
x = 0.33, 0.4, 0.5, 0.6 and 0.75). 

The detailed structures were fitted via Rietveld refine- 
ment using the FULLPROF program [28 . Fig. 3 (upper) 
presents the refined neutron powder diffraction data col- 
lected at HB-2A for the x = 0.1 sample at the base tem- 
perature of 4K. Fig. 3 (lower) shows a closer look at 
the evolution of the 220 nuclear peak for x = 0.1, 0.33 
and 0.75. The structural transition is evident in the x 
= 0.1 sample as the peak splitting in 220at is easily dis- 
cernible; this is much less obvious for samples with x > 
0.33. Fig. 4 shows the temperature dependence of the 
relative goodness of fits obtained by fitting the full pow- 
der patterns to each of the two different known structural 
models (tetragonal P4/nmm vs. orthorhombic Cmma.) 
It is apparent that the orthorhombic model fits better 
than the tetragonal model for x = 0.1 at T < 120K and 
for x = 0.33 at T < 75K. The inferred structural phase 
transition temperature in the x = 0.33 sample is consis- 
tent with the anomaly near 80K observed in the resis- 
tivity datapj)]. The diffraction measurements for higher 
Ru concentrations provide no evidence for any structural 
transition, suggesting that it is suppressed for Ru doping 
between 33-40%. 

Interestingly for other isoelectronically doped 1111 ma- 
terials the structural phase transition has been observed 
to persist to dopings well above the 10% level; for ex- 
ample in CeFeAsi-^P^O both the transition to an or- 
thorhombic structure and antiferromagnetic order are 
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TABLE I. Structural parameters and x 2 values from Rietveld refinement of NPD data from HB2A (for Ru% = 10 only) 



and POWGEN, at T = 12K. (Space 


group: Cmma (#67) for 10% and 33%, P4/nmm (#129) otherwise.) 




Ru% a/b (A) 


c(A) 




Zp r 




x 2 


10 3.96865(8)/3.98509(8) 


8.5622(2) 


0.6550(5) 


0.1402(7) 


71.89(9)/72.23(9) 


3.34 


33 4.00063(3)/4.00733(3) 


8.5154(1) 


0.6570(3) 


0.1389(5) 


71.95(7)/72.09(7) 


3.64 


40 4.01162(4) 


8.4911(2) 


0.6582(5) 


0.1378(6) 


71.97(7) 


3.51 


50 4.02491(6) 


8.4613(2) 


0.6626(2) 


0.1328(5) 


71.43(7) 


2.82 


60 4.03351(7) 


8.4364(3) 


0.6603(6) 


0.1333(9) 


71.93(6) 


4.04 


75 4.05088(4) 


8.3980(1) 


0.6587(4) 


0.1366(6) 


72.41(7) 


3.26 
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FIG. 5. (a): The intensity of the 101m peak as a function 
of temperature. The AFM ordering of both Pr and Fe are 
visible in x = 0.05 but not x = 0.1. The intensity of the 
101m peak reaches a domelike maximum around 60K for x 
= 0.05. Solid dots represent single data points collected from 
counting at the center peak position; hollow dots represent 
normalized data from fits of entire magnetic peaks. The lines 
and shadings are guides to the eye; (b): Intensity of the 100m 
peak as a function of temperature, showing the size difference 
of the Pr moment between x = 0.05 and 0.1. 
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FIG. 6. Neutron powder diffraction data collected from POW- 
GEN refined to the orthorhombic Cmma structure for x = 
0.33 at T = 12K. 



present up to x = 0.4[16 . In contrast, a much smaller 
level of carrier doping is needed to affect the transitions. 
In RFeAsOi-^F^ all structural and magnetic transitions 
were suppressed with less than 10% doping [l] l29ti32] . In 
fact the doping level resulting in the highest observed su- 
perconducting transition temperature occurs at x = 0.1 
for R = La[T], 0.11 for Pr[29l|30]. 0.11 for Nd[31] and 0.1 
for Sm[32], with superconductivity usually disappearing 
around x = 0.2 [30H32]. 

We now turn to a detailed examination of the magnetic 
transition that is observed in PrFe i-^Ru^ AsO with 5% 
Ru doping. In Fig. 5 we present temperature depen- 
dent measurements of the intensities of the 101m and 
100m magnetic peaks measured using HB-lA for doping 
levels x = 0.05 and 0.1. The 101m peak shows contri- 
butions from both Pr and Fe moments while the 100m 
peak is sensitive only to the Pr moments. It is apparent 
in Fig. 5(a) and (b) that magnetic ordering of both Pr 
and Fe exists at x = 0.05 but is suppressed at x = 0.1. 
Fe orders antiferromagnetically in the x = 0.05 sample 
at 120K which may be compared to 137K[27 observed 
in the undoped PrFe AsO; no AFM order is visible in x = 
0.1. Evidently in PrFe i-^Ru^ AsO, the magnetic transi- 
tion is more sensitive to Ru doping than is the structural 
transition. In any case neither transition is observed for 
x > 0.4. 

The temperature dependence of the magnetic inten- 
sity of the 101m peak is particularly interesting. As seen 
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in Fig 5(a), in the x = 0.05 sample the intensity of the 
101m peak reaches a dome-like maximum at T ~ 60 K, 
followed by a modest decrease as the temperature is low- 
ered, increasing again below T ~ 14 K, the AFM ordering 
temperature for Pr. This feature has also been observed 
in undoped PrFeAsO at T - 60K[24 . Critical fluctu- 
ations in the Pr magnetic subsystem may contribute to 
this phenomenon [24] . 

To examine the temperature dependence of the crys- 
tal structure and lattice constants and in particular 
the behavior of the NTE we performed higher res- 
olution neutron powder diffraction at POWGEN for 
PrFei-^Ru^AsO with x = 0.33, 0.4, 0.5, 0.6 and 0.75. 

Fig. 6 illustrates a typical refinement using the pow- 
der pattern for x = 0.33. Table I shows a summary of 
the lattice and crystallographic parameters obtained at 
12 K for several different concentrations of Ru. The an- 
gle (j) connecting near-neighbor Fe atoms via an As atom 
(see Fig. 1) is also included. To within the resolution 
of the experiment for each concentration the crystallo- 
graphic parameters zas, zp ri and the angle <p have only 
a very small temperature dependence below 200 K. It 
should be noted however that the precision with which 
these parameters are determined is an order of magni- 
tude coarser than the precision with which the lattice 
constants are determined. This is expected since the lat- 
tice constants are determined precisely by the peak po- 
sitions alone whereas the other parameters are refined 
by fitting peak intensities. On the other hand there is a 
clearly observable dependence of the lattice constants on 
the temperature over the entire measured range. 

Fig. 7 (upper) shows the contrast in the behavior of 
the a- and c-axes' doping dependence at 12K. The trends 
are consistent with the room temperature x-ray data re- 
ported by McGuire et at. [19], which they noted may be 
related to both the larger radius of Ru compared to Fe as 
well as constraints on Ru-As and Fe-As bond distances. 
Fig. 7 (lower) shows the temperature dependence of the 
lattice parameters for the concentration x = 0.4. The 
dashed red line is the extrapolated behavior for the c- 
axis based on a Gruneisen fit to the volume defined by c 3 
for the data at 40K and above. This is a representative 
concentration as the others show similar behavior. At 
low temperatures (i.e. for T < 50K), the lattice param- 
eters deviate from the Gruneisen model fit in opposite 
directions, somewhat conserving the overall cell volume. 
This is illustrated in Fig. 8, which shows the temperature 
dependence of the unit cell volume. A Gruneisen model 
fits the entire cell volume better than it fits either the a- 
or c-axis alone. This behavior is not uncommon among 
materials that are anisotropic, where expansion along one 
direction is often accompanied by contraction along an- 
other in order to preserve the overall volume; graphite is 
one well known example [33 . Despite this compensation 
a modest NTE of the unit cell volume is still visible in 
all samples for T < 20K. 

NTE in the c-axis is a prominent feature observed 
previously in the undoped parent compound: the c-axis 
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FIG. 7. (Upper): Doping dependence of a and c at 12 K; 
(Lower): Temperature dependence of a- and c-axis plotted 
together, showing the contrasting thermal behavior between 
the two. The red line is the Gruneisen model prediction for 
the c axis based on fits to the data for T > 40K as described 
in the text. 
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reaches a minimum at T = 50K, and then expands as the 
system is further cooled [24]. Kimber et al. ascribed the 
presence of c-axis NTE in the non-superconducting par- 
ent compound to the effects of spin-lattice coupling [24J. 
This was based on their observation that both mag- 
netic order and NTE are absent in the superconduct- 
ing PrFeAsOi-^F^, and the proximity of the tempera- 
ture dependence of the NTE and that of the dome-like 
feature observed in the intensity of magnetic peaks as- 
sociated with the ordering of the Fe moments [24j. In 
this context, the observation of NTE in PrFei-^Ru^ AsO 
with x > 0.1 is surprising, since long range magnetic or- 
der is absent. Fig. 9 shows the temperature dependence 
of the fractional change in the refined lattice parameters 
a, b and c for several concentrations. As seen clearly in 
Fig. 9 (upper), the anomalous NTE along the c-axis is 
observed in all samples up to at least x = 0.75. Figure 
10 shows a phase diagram of PrFei-^Ru^AsO with the 
region of c-axis NTE indicated. 

In contrast the thermal behavior of the a-axis is almost 
linear at low temperatures, with the fractional change al- 
most independent of the Ru concentration. The existence 
of c-axis NTE in concentrations without magnetic order 
suggests that the origin of this behavior might involve 
more than simply spin-lattice coupling. Although there 
is little published data on the c-axis temperature depen- 
dence of superconducting 1111 materials, the NTE is ab- 
sent among those for which data is accessible, including 
PrFeAsO0.s5F0.15 [24 1, and LaFeAsOi_^F^p4l [35]. To 
date we know of no example in the literature of a su- 
perconducting 1111 compound that does show NTE. 

It is interesting to compare the low temperature ther- 
mal variation of the a- and c-axes with the composition 
dependence shown in Fig. 7. In both cases the param- 
eters vary in the opposite direction. The reason for the 
opposite temperature dependence of the a- and c-axes 
remains an open question. The combination of NTE in 
the c axis and near conservation of volume may lead to 
subtle but important changes in the effective magnetic 
and electronic interactions of the Fe-As layer both inter- 
nally and with the rest of the sample. Therefore, it can 
possibly be inferred that the NTE is linked to the lack of 
superconductivity (or any new ground state) in the Ru 
doped samples. 



FIG. 9. Detailed temperature dependence of of the lattice 
constants for several values of Ru doping. (Upper): The quan- 
tity Ac/cmok is defined as (c(T)-c(T=140K))/c(T=140K). 
NTE is clearly visible in all samples for T < 50K. The solid 
line is a representative Gruneisen model fit for x = 0.33 also 
shown in Fig. 8. (Middle): A similar plot, showing the rela- 
tive change in a- and b-axis with respect to lattice parameter 
a measured at 140K, which is above the structural transition 
for x = 0.1 and 0.33; (Lower): A similar plot for a-axis for x 
= 0.4, 0.5 and 0.75, showing the lack of NTE at low temper- 
atures. 



IV. CONCLUSION 

We have found that magnetic order on both Fe and 
Pr sites in Pr-1111 is suppressed with 10% of isoelec- 
tronic doping by Ru/Fe substitution, and the tetragonal- 
orthorhombic structural transition with 40%. We have 
studied the effect of Ru doping on the anomalous NTE 
along the c-axis and found that the c-axis NTE is present 
in all of our samples. For all concentrations the lattice 
parameter along the c-axis reaches a minimum around 
50K. The NTE is extremely resilient to Ru doping and 
is measurable up to at least x = 0.75. The c-axis NTE 
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Ru concentration (x) 

FIG. 10. Phase diagram of PrFei_ x Ru x AsO constructed with 
data reported in this article, showing temperatures for the 
tetragonal-orthorhombic transition (T structural), the mag- 
netic Fe ordering (Tsdw), and the magnetic Pr ordering 
(Tp r ), in the measured Ru concentrations. The region of 
c-axis NTE presents in all Ru concentrations below T = 50K. 



extends across the entire composition range studied: in 
compositions which remain tetragonal to the lowest tem- 
peratures, those which exhibit tetragonal-orthorhombic 
distortions with no long-range magnetic ordering, and 



those that undergo both structural and magnetic phase 
transitions. Other anomalous behaviors observed include 
the near linear temperature expansion of a-axis at low 
temperature, which also has little discernible effect from 
doping. The contrasting thermal behavior between a- 
and c-axis somewhat conserves the unit cell volume. The 
temperature dependence of the unit cell volume agrees 
better with the Gruneisen prediction than the a- or c- 
axis alone. It would be interesting to investigate whether 
the apparent relationship of NTE to the lack of super- 
conductivity is indeed a general rule for 1111 materials. 
Such knowledge should provide important insights into 
the underlying mechanism for superconductivity in Fe 
based materials. 
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